Известия Томского политехнического университета. 2006. Т. 309. № 8 


УДК 535:621.372.8 


ДИНАМИЧЕСКИЕ ИСКАЖЕНИЯ СИГНАЛА В КОРОТКИХ ОТРЕЗКАХ 
МНОГОМОДОВЫХ ОПТИЧЕСКИХ ВОЛОКОН 

В.Э. Аппельт, А.С. Задорин, Р.С. Круглов, П.В. Чернов 

Томский государственный университет систем управления и радиоэлектроники 
Е-таіІ: Апа{оІу.2асІогіп@г2І.{и$иг.ги 

Работа посвящена моделированию процессов линейных динамических искажений оптических сигналов в многомодовых волок- 
нах, обусловленных дисперсией волноводных мод и межмодовой связью, вызываемой шероховатостями поверхности сердце- 
вины оптического волокна. Показано, что при длинах оптического волокна, превышающих длину установившейся вязи мод, сиг- 
нал приобретает устойчивую форму, не зависящую от условий возбуждения волокна, и близкую к отклику электронной интегри- 
рующей цепи. Проведено сравнение результатов расчетного и натурного экспериментов по исследованию зависимости ширины 
полосы пропускания от длины волокна в полимерном оптическом волокне со ступенчатым профилем показателя преломления, 
демонстрирующее хорошее соответствие расчетных и экспериментальных данных. 


Введение 

Многомодовые оптические волокна (МОВ), 
как правило, используются при построении систем 
связи, радиус которых настолько мал, что длины 
элементарных кабельных участков соизмеримы с 
длиной установившейся связи мод. Математиче- 
ские модели, описывающие трансформацию опти- 
ческих импульсов в указанных условиях, развиты 
слабо, что накладывает ограничение на возможно- 
сти оптимизации быстродействия телекоммуника- 
ционных сетей. 

В оптически линейных МОВ основным меха- 
низмом искажений сигнала является эффекты 
межмодовой дисперсии (ММД) [1, 2], а трудности 
их моделирования связаны, прежде всего, с нали- 
чием микронерегулярностей на границе раздела в 
среднем ровной поверхности сердцевины оптиче- 
ского волокна (ОВ). Такие шероховатости сердце- 
вины приводит к формированию многократных 
связей внутри спектра волноводных мод (ВМ), и 
радиационным потерям. 

Основной целью настоящей работы является 
моделирование процессов трансформации оптиче- 
ских сигналов в МОВ, обусловленных дисперсией 
волноводных мод и межмодовой связью, вызы- 
ваемой шероховатостями поверхности оптическо- 
го волокна. 


Постановка задачи 


Рассмотрим слабонаправляющее, И- модовое, 
оптически линейное ОВ [3] со ступенчатым про- 
филем показателя преломления, без потерь, в на- 
правлении оси і которого распространяются ли- 
нейно-поляризованная ВМ (ЬР) [1]: 
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где N 1 - функции Бесселя и Неймана порядка /, 
описывающие профиль ВМ в области сердцевины 


и оболочки соответственно; е, - вектор поляриза- 
ции ВМ; Е, - амплитуда ВМ; со - круговая частота; 
Д - постоянная распространения ВМ; а - коэф- 
фициент затухания; г, ф, і - координаты цилин- 
дрической системы координат; 

« = а ^ к 1 п 1о ~ Я 2 > и; = - кЛ , 

где а - радиус сердцевины, п со и п с , - показатели 
преломления сердцевины и оболочки соответ- 
ственно, А 0 - волновое число световой волны в ва- 
кууме. 

Среднюю интенсивность всех мод обозначим 
вектором интенсивностей I. Компоненты I будем 
называть родительскими волноводными модами. 
Будем считать, что вектор I модулирован динами- 
ческой зависимостью х(1). Обозначим соответ- 
ствующий модовый пакет х(1)- 1 вектором х(7). 
Входной сигнал оптической системы передачи х(1), 
распространяясь вдоль оси 1 0В, под действием ря- 
да физических механизмов претерпевает линейные 
искажения, трансформируясь к виду 

Решение задачи моделирования и снижения 
линейных ММД-искажений сигнала в регулярном 
ОВ с идеальной поверхностью сердцевины не вы- 
зывает затруднений, т. к. в указанной ситуации 
межмодовая связь исключена, и интерференция 
ВМ возможна только на фотоприемнике приемно- 
го модуля. 

Рассмотрим далее влияние на у (і,1) разного ро- 
да отклонений от указанной идеализации. Среди 
таковых ниже рассматриваются микроскопические 
неоднородности ОВ, сосредоточенные на поверх- 
ности сердцевины. 

Матричная импульсная характеристика МОВ 

Микроскопические неоднородности в среднем 
ровной поверхности сердцевины могут быть опи- 
саны изотропным двумерным стохастическим по- 
лем Дг). Радиус корреляции поля обозначим через 
4, а среднюю высоту шероховатости - 8. 

Для спектра І 0 распространяющихся в напра- 
влении оси г волноводных мод локальные микро- 
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скопические неоднородности поверхности сердце- 
вины эквивалентны пространственным флуктуа- 
циям Аи эффективного показателя преломления 
п(г) относительно его среднего уровня п ср [4]. 

Указанные шероховатости О В приводят к ра- 
спределенной по всему волокну связи внутри спек- 
тра ВМ и непрерывному и многократному энерго- 
обмену между составляющими сигнального векто- 
ра у при его распространении вдоль оси і. 

При анализе параметров системы связи практи- 
ческий интерес представляет интенсивность цифро- 
вого сигнала, усредненная по его длительности г. 
Для реальных систем с многомодовыми ОВ это вре- 
мя значительно превышает время прохождения сиг- 
налом участка длиной І с . Поэтому, в соответствии с 
[5], мы будем рассматривать лишь накапливающие- 
ся изменения интенсивностей у усредненные 
по интервалу А, значительно превышающему /„ но 
существенно меньшему, чем длина нормализации Ь к 
к - й «родительской» волноводной моды (РВМ), на 
которой спектр І(Д трансформируется в устойчи- 
вый, равновесный набор мод, т. е. 

1 С « А«ѵт«Ь ь (2) 

где ѵ - групповая скорость волнового пакета. 

Отклик у(Д,/) на выходе элементарного участка А 
может быть определен как свертка импульсной ха- 
рактеристики с входным сигналом х(г). Однако в слу- 
чае многомодового волокна импульсная характери- 
стика Ь(?) является матричной величиной. Далее бу- 
дем назвать ее импульсной матрицей межмодового 
рассеяния (ИММР) [6]. По своему значению ИММР 
аналогична многомерной импульсной характеристи- 
ке линейного электронного многополюсника [7]. 

Для расчета динамики у (?) на выходе участка ОВ 
с длиной кратной А, т. е. ^ = рА, воспользуемся 
свойством оптической линейности ОВ. Последова- 
тельно р раз применим свертку импульсной харак- 
теристики с входным сигналом х(?) 

у(/>-А,0 = Ь(0*...* х(0, (3) 

Р 

где символом *...* обозначена /ь кратная матричная 
временная свертка сигнала х(/) с ИММР. 

Очевидно, что размерность й(/) определяется 
числом ВМ, в которые рассеивается энергия векто- 
ра РВМ х(/). Если принимать во внимание только 
те моды, которые вследствие рассеяния энергии 
РВМ, распространяются в попутном направлении, 
то размерность указанной матричной импульсной 
характеристики оказывается равной ЫхК Наличи- 
ем обратно рассеянных мод можно пренебречь, т. к. 
их влияние сказывается только лишь на ближнем 
конце системы передачи. 

Импульсная характеристика любого линейного 
многополюсника, представляет собой динамиче- 
ский отклик системы на сигнал, близкий к (5-функ- 
ции. В работах [8, 9] было установлено, что элемен- 
ты импульсной матрицы межмодового рассеяния 
й(/) определяются выражением 


к тіп = (М што / т т )Кесі((1 - т т )/т т ) ~ М шт 8(1 - т т ), при 
к тп = (М тп / т тп )Кесі((і - т п ) / т тп ) , при 
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где т ш - относительная временная задержка между 
я-й ВМ «потомком» (ПВМ) и т - й РВМ, вычисля- 
емая через разность групповых скоростей ѵ„ и ѵ т ука- 
занных мод; М - матрица межмодового рассеяния 
(ММР), определяяемая следующим выражением [8] 
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где О тп (Ак) - двумерный дискретный энергетиче- 
ский спектр шероховатости; - коэффициент 
связи между т - й и я-й модами: 

оо 

$.„ = /а 8 г (г)-Е:(г).Е:(г)Ыг, (7) 


Е' - профиль ВМ; Ае=я-Ая - вариация диэлектри- 
ческой проницаемости среды в поле простран- 
ственной гармоники А к=к,-к т энергетического 
спектра шероховатостей С т „(Ак), обеспечивающей 
синхронизм между волноводными модами. 

Коэффициент затухания а„ входящий в (5, 6), 
определяет радиационные потери ВМ, обусловлен- 
ные наличием шероховатости поверхности сердце- 
вины [7]. 

Соотношение (3), вместе с (4-6) формализуют 
эффект ММД и имеют прозрачный физический 
смысл. На элементарном участке А многомодового 
ОВ любая из N родительских мод испытывает за- 
держку т т =А/ѵ т и формирует /V— 1 мод потомков. 
При этом, в силу межмодовой задержки, перекачка 
энергии от РВМ к ПВМ оказывается рассредото- 
ченной во временном интервале т ір что, в полном 
соответствии с (4), указывает на динамическое 
уширение формы ПВМ относительно (5-импульса 
РВМ. ПВМ описываются прямоугольными им- 
пульсами, которые могут как отставать, так и опере- 
жать родительский 5-импульс. Представленная мо- 
дель описывает ММД-искажения оптического сиг- 
нала лишь при малых длинах ОВ А и слабом рассе- 
янии энергии РВМ. При больших длинах ОВ тран- 
сформация импульса описываются многократной 
матричной сверткой (3) вектора РВМ с ИММР, ко- 
торая является физической моделью эффекта мно- 
гократного рассеяния и интенсивного перемешива- 
ния мод в ОВ. В ходе этого процесса вектор РВМ 
х(/) испытывает как радиационное затухание, так и 
сильное взаимодействие с ПВМ [8, 9]. 


Результаты численного моделирования и их анализ 

Численное моделирование динамики тран- 
сформации формы импульса проводилось на осно- 
ве выражений (3), а также формул (4-6), которые 
позволяют определить все элементы ИММР. 
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Основные параметры ОВ: показатели прело- 
мления сердцевины и оболочки « го =1,504, 
и с ,=1,500; радиус сердцевины а= 9,7 мкм, при этом 
на длине волны излучения А=650 нм обеспечива- 
лось распространение 16 ВМ (7Ѵ=15). Параметры 
шероховатости поверхности сердцевины: радиус 
корреляции 1=5 мкм; средняя высота шероховато- 
сти 5=5 нм. 

Результаты моделирования трансформации 
формы импульса при возбуждении ОВ мода- 
ми различного порядка представлены на рис. 1,2. 

Форсированное искажение формы переднего 
фронта при возбуждении ОВ модой ЬР 6І обусловле- 
но в первую очередь перекачкой энергии из роди- 
тельской ВМ в ПВМ более низкого порядка, груп- 
повая скорость которых выше. После того, как 
спектр ВМ в процессе самофильтрации достигает 
равновесного состояния, форма импульса у(1,р) 
становиться близкой к отклику электронной инте- 
грирующей цепи. 


При возбуждении ОВ модой более низкого по- 
рядка ЬР п энергия РВМ на начальных участках ОВ 
эффективно передается к ПВМ более высокого по- 
рядка, что проявляется в затягивании заднего фрон- 
та импульса. На расстояниях, превышающих длину 
нормализации модового спектра, форма импульса 
у(1,і) не зависит от условий возбуждения ОВ, т. е. 
можно говорить о нормализации формы импульса. 

Результат моделирования среднеквадратичного 
уширения [10] импульсной характеристики ОВ при 
его возбуждении модой ЬР т представлен на рис. 3. 

Среднеквадратичное уширение импульсной ха- 
рактеристики при возбуждении основной модой 
линейного многомодового ОВ на малых длинах 
пропорционально дайне ОВ Ь, а при больших - 
пропорционально УІ. 

В литературе, как правило, рассматриваются 
лишь предельные случаи поведения представленной 
кривой [10]. Если длина ОВ превышает длину рав- 
новесного модового распределения, то спектр вол- 



і/т 

Рис. 1. Нормализация оптического сигнала с прямоугольной огибающей у( і,г) в ОВ при его возбуждении модой ІР е 1 (к=11) 



1/т 

Рис. 2. Нормализация оптического сигнала с прямоугольной огибающей у( і,г) в О В при его возбуждении модой /Р 31 (к=4) 
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новодных мод сохраняет свое равновесное состоя- 
ние (состояние нормализации), при этом сг~Ѵі. В 
том случае, если длина ОВ значительно меньше дли- 
ны равновесного модового распределения, то влия- 
ние межмодовой связи на зависимость аЩ оказы- 
вается пренебрежимо малым, поэтому используется 
следующая асимптотическая оценка: а~Ь. 



I I I I I 11111 I I I I ■ ...I I I I I I і Гі ■ ■ ■ ■ 1 і ГІ1І 

10 100 1000 
Рис. 3. Дисперсионное уширение сигнала о(і) в МОВ при 
возбуждении моды ТР 0 1 


При моделировании динамических искажений 
сигналов, распространяющихся по реальным мно- 
гомодовым оптическим волокнам, остро встает 
проблема нехватки вычислительных ресурсов. В 
полимерном оптическом волокне (ПОВ), радиус 
сердцевины которого й=490 мкм, числовая аперту- 
ра Жу 4=0,46, на длине волны А=650 нм возможно 
распространение более чем 2,37' ІО 6 ВМ. При этом 
ИММР для указанного ОВ будет состоять пример- 
но из 5,63- ІО 12 элементов, каждый из которых в 
свою очередь является вектором, описывающим 
импульсную характеристику соответствующей ВМ. 

Задача аппроксимации спектра ВМ в первую 
очередь связана с проблемой выбора профиля фор- 
мальной волноводной моды, которая должна заме- 
нить собой группу реальных ВМ. Амплитуда ВМ в 
области шероховатости определяет уровень радиа- 
ционных потерь и коэффициент связи между ВМ. 

В силу того, что для ПОВ характерен относи- 
тельно большой радиус сердцевины, амплитуды 
пространственных угловых (азимутальных) гармо- 
ник энергетического спектра шероховатости прак- 
тически равны. Поэтому эффективность взаимо- 
действия различных ВМ, а также величина радиа- 
ционных потерь ВМ будут определяться только ра- 
диальным порядком моды. 

Расчет параметров формальных волноводных 
мод удобно проводить с помощью характеристиче- 
ского уравнения, которое устанавливает связь меж- 
ду азимутальным / и радиальным т порядками ВМ 
и ее постоянной распространения Д [ 1 1 ] 

«К 2 Л 2 -^-ф 2 ] Х - 


-/ С08 - 


<^ 2 Л 2 -Д 2 -ф 2 



где т - радиальный порядок, указывающий коли- 
чество узлов поля ВМ в радиальном направлении; 
ф а - фазовый сдвиг, возникающий при полном вну- 
треннем отражении от стенки сердцевины ОВ. 

При формировании ИММР, описывающей ди- 
намику трансформации импульса, переносимого 
группой формальных волноводных мод, также 
необходимо ввести дополнительный дисперсион- 
ный коэффициент уширения т' т для каждой роди- 
тельской формальной волноводной моды. Указан- 
ный коэффициент позволит учесть эффекты меж- 
модовой дисперсии, возникающие внутри каждой 
аппроксимированной группы реальных ВМ. 

Дисперсионный коэффициент уширения роди- 
тельских формальных волноводных мод т’„, можно 
записать в виде 

( 8 ) 

где Ѵ т и Ѵ т - групповые скорости реальных ВМ, 
находящихся на противоположных краях фрагмен- 
та спектра ВМ, аппроксимируемого одной фор- 
мальной волноводной модой. 

С учетом (8) выражение, определяющее матри- 
цу ИММР (4), может быть переписано в виде 

К т =( м ™/<Жес1;((Т-т т )/<), при »;=«,] 

Кп = ( м „,„ /г т „)ВссІ((/-г„)/г„ ш ) , при тФѣу 

где Кесі - функция, описывающая импульс с пря- 
моугольной огибающей и длительностью г. 

Зависимость ширины полосы пропускания 
ПОВ от длины волокна, рассчитанная с примене- 
нием описанного выше алгоритма аппроксимации 
спектра ВМ, представлена на рис. 4. 

Основные параметры ПОВ, использованные в 
численном моделировании, соответствуют стан- 
дартному ПОВ Тоту РТИ-СОКЮТ. показатели пре- 
ломления сердцевины и оболочки п с = 1,492, 
п с і =1,419; радиус сердцевины а=490 мкм. Аппрок- 
симация спектра ВМ осуществлена 16-ю формаль- 
ными волноводными модами. Длина волны излу- 
чения А=650 нм. 

Кружками на рис. 4 показанные эксперимен- 
тальные данные для ПОВ Тоту РТИ-С01001. 

При моделировании зависимости, представлен- 
ной на рис. 4, параметры шероховатости поверхности 
сердцевины были выбраны таким образом, чтобы 
обеспечить минимальное расхождение между экспе- 
риментальными и расчетными данными. Получен- 
ные значения параметров шероховатости поверхно- 
сти сердцевины ОВ (радиус корреляции 1=2 мкм; 
средняя высота шероховатости <5=300 нм). Незначи- 
тельное расхождение полученных результатов об- 
условлено наличием дополнительных, неучтенных в 
модели, факторов, приводя щ их к межмодовому взаи- 
модействию и дополнительным радиационным поте- 
рям (например, рассеяние Рэлея). 

На начальных участках, значительно меньших 
длины установившейся вязи мод, и на длинах, зна- 
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чительно превышающих эту характерную дистан- 
цию, поведение зависимости ширины полосы про- 
пускания от длины О В соответствует хорошо из- 
вестным асимптотам [10]. 



10 100 1000 10000 


Рис. 4. Зависимость ширины полосы пропускания от длины 

пов 

Особенностью каждой из представленных на 
рисунке 4 ветвей зависимости АД V) является нали- 
чие на них участков, на которых значение про- 
странственной производной сІА/Щ/сИ достигает 
своего максимума, превышающего ее уровень в 
асимптотических областях. В пределах указанных 
участков происходит наиболее быстрое, форсиро- 
ванное снижение полосы пропускания О В при уве- 
личении его длины. Объяснение данного эффекта 
не вызывает затруднений. Действительно, именно 
в указанных областях процесс нормализации модо- 
вою спектра ОВ переходит в свою завершающую 
фазу. Первоначально наиболее мощные рабочие 
моды ослабляются и рассеиваются в ПВМ на шеро- 
ховатостях предшествующих участков ОВ. В ука- 
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занных областях они сравниваются по мощности с 
остальными «модами-потомками». Естественно, 
что в условиях переноса сигнала модовым спек- 
тром с равномерным распределением амплитуд, 
уровень, обусловленных межмодовой дисперсией, 
линейных искажений сигнала оказывается макси- 
мальным. 

Заключение 

Представлена математическая модель, описы- 
вающая эффекты линейных искажений сигнала в 
многомодовых оптических волокнах при наличии 
неоднородностей поверхности сердцевины. 

Проведено сравнение результатов расчетного и 
натурного экспериментов по исследованию зави- 
симости ширины полосы пропускания полимерно- 
го оптического волокна со ступенчатым профилем 
показателя преломления от его длины, демонстри- 
рующее хорошее соответствие расчетных и экспе- 
риментальных данных. 

На длинах оптического волокна, значительно 
меньших длины установившейся связи мод, либо 
значительно превышающих ее, поведение зависи- 
мости ширины полосы пропускания волокна от его 
длины соответствует хорошо известным асимпто- 
там. 

Установлено, что отклик нерегулярного много- 
модового оптического волокна на расстояниях, 
превышающих длину установившейся вязи мод, 
стремится к форме отклика электронной интегри- 
рующей цепи независимо от условий возбуждения 
оптического волокна. 

Работа выполнена при финансовой поддержке организа- 
ции ШТА8 (грант № 04-83-3239). 
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